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We introduce a multitype contact process with temporal hetero- 
geneity involving two species competing for space on the d-dimensional 
integer lattice. Time is divided into seasons called alternately season 
1 and season 2. We prove that there is an open set of the parame- 
ters for which both species can coexist when their dispersal range is 
large enough. Numerical simulations also suggest that three species 
can coexist in the presence of two seasons. This contrasts with the 
long-term behavior of the time-homogeneous multitype contact pro- 
cess for which the species with the higher birth rate outcompetes the 
other species when the death rates are equal. 

1. Introduction. Consider an ecosystem inhabited by multiple species 
competing for a single resource represented by a function R. Let ui, i = 
1,2, ... ,K, denote the population density of species i. In the original models 
of Lotka (1925) and Volterra (1928), the per capita growth rate for each 
species increases linearly with the amount of resource available: 

(1) — < ^=(3iR(u 1 ,u 2 ,...,u K )-5i, 
Ui at 

where Pi and 5i, respectively, denote the birth and death rates of species 
i. When R is a decreasing function of population densities, Volterra (1928) 
showed that only the species with the highest Pi to 5i ratio survives as 
t — > oo. 

When R = 1 — {u\ + + • • • + uk) denotes the density of space available, 
the Lotka-Volterra model (1) can be seen as the mean-field approximation, 
or nonspatial deterministic analogue, of the multitype contact process in- 
troduced by Neuhauser (1992). Assuming for simplicity that the number of 
species is K = 2, her model is a continuous-time Markov process whose state 
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at time t is a function r]t that maps the d-dimensional integer lattice Z into 
{0,1,2}, with denoting empty sites, and 1 and 2 denoting sites occupied 
by an individual of species 1 and 2, respectively. The evolution at site x G Z d 
is described by the following transition rates: 

at rate /3i p(y,%) ^-{n(y)=i\i 1 — > at rate<5i, 
at rate 37)1/^(^=2} , 2 — > at rate 62, 

where p(y,x) is a translation invariant transition probability on Z rf . This in- 
dicates that species of type i, i = l,2, produces offspring at rate and dies 
at rate 5{. Moreover, an offspring produced at site y G Z d is sent to a site x 
chosen at random with probability p(y, x). If site x is vacant, it becomes oc- 
cupied by the offspring, otherwise the birth is suppressed. Neuhauser (1992) 
proved for the multitype contact process that if <5i = 62 then the species 
with the larger birth rate outcompetes the other species. She also conjec- 
tured that, as predicted by the model (1), the species with the larger /3$ to 
6i ratio takes over. While it is believed that there is no coexistence in the 
Neuhauser's multitype contact process, there are a number of interacting 
particle systems for which coexistence has been proved rigorously. 

One mechanism that allows coexistence is when intraspecific competition 
exceeds interspecific competition. This has been proved by Neuhauser and 
Pacala (1999) for a voter type model in which each site of the lattice is 
occupied by an individual of type 1 or 2. When both species have the same 
birth rate, the state at site x G Z rf evolves according to the transition rates 

l->2 at rate/ 2 (/i + CK12/2), 
2^1 at rate /i(/ 2 + a 2 i/i), 

where /j = f%{x) denotes the fraction of neighbors of site x in state i. The pa- 
rameters CKjj model the pressure of species i on species j. Except for the 
1-dimensional nearest neighbor case, in the absence of interspecific competi- 
tion, that is, a 12 = 0-21 = 0, there is a stationary distribution in which sites 
are independent [see Theorem 1 in Neuhauser and Pacala (1999)], indicating 
that both species coexist. While increasing a\2 and CK21, spatial correlations 
build up, but coexistence is maintained as long as the interspecific com- 
petition is sufficiently small. The inclusion of a spatial structure, however, 
translates into a reduction of the set of parameters for which coexistence 
occurs. 

A second mechanism that promotes coexistence is environmental hetero- 
geneity. Suppose that the regular lattice Z 2 is inhabited by two types of 
plants. The environment is static and described by a partition {Hi, #2} of 
the space with Hj denoting the set of sites occupied by plant i. Consider two 
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specialist symbionts; for i = 1,2, symbiont of type i feeds on plant i. This 
can be modeled by a multitype contact process in which births of symbionts 
of type i are suppressed whenever the offspring is sent to a site x € Hj with 
i ^ j. Since the habitats H± and H2 are disjoint and so symbionts 1 and 2 do 
not interact, the problem of coexistence reduces to the following percolation 
problem: If for i = 1, 2, there exists an infinite sequence of sites in Hi along 
which symbiont i can spread then coexistence is possible. This depends on 
both the spatial arrangement of Hi and the dispersal range of symbiont i. 

Lanchier and Neuhauser (2006) proved that coexistence occurs in a non- 
trivial case in which the habitats of both symbionts overlap. Their model 
is referred to as the static-host model. Symbiont 1 is a specialist feeding on 
plant 1 while symbiont 2 is a generalist feeding on both plants. The universe 
is tiled with N x N squares. Whether a square is occupied by plants of type 

1 or 2 is determined by flipping independent fair coins. In this situation, 
coexistence occurs when the birth rate of the specialist is larger than that 
of the generalist and the square size N is much larger than the dispersal 
range of the symbionts. In the mean- field approximation of the static- host 
model, the birth rate of the specialist needs to exceed twice the birth rate 
of the generalist so that coexistence occurs. Lanchier and Neuhauser (2006) 
also introduced a stochastic process in which spatial heterogeneities are no 
longer static but generated by the symbionts. This model is referred to as 
the dynamic-host model. In the absence of symbionts, plants of type 1 and 

2 compete according to a biased voter model with a selective advantage for 
the l's. In particular, plants of type 1 outcompete plants of type 2. Intro- 
ducing a specialist pathogen feeding on plants of type 1 whose harmful effect 
is modeled by an increase of the death rate of infected plants, Durrett and 
Lanchier (2008) proved that coexistence may occur provided the dispersal 
range of the species is sufficiently large, but conjectured that coexistence 
is not possible if the symbiont is a specialist mutualist feeding on plants of 
type 2. 

In some cases, the presence of two types of space is not immediately obvi- 
ous. Chan and Durrett (2006) considered a generalization of the multitype 
contact process in which species 1 is a good competitor with a short range 
dispersal kernel, whereas species 2 has a long range dispersal kernel which 
makes it a good colonizer. In this situation, occurrences of forest fires (i.e., 
removal of all the individuals contained in a large square) at an appropriate 
rate allows species 2 to survive by migrating to newly created patches. In 
this model, there are two types of space: blocks freshly cleared by a fire 
which are occupied by individuals of species 2 and older blocks which are 
occupied by individuals of species 1. 

The aim of this article is to prove that temporal heterogeneity can also 
promote coexistence in a stochastic spatial model. Armstrong and McGehee 
(1976) have shown in an ODE model that coexistence can occur for species 
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living in a periodic environment with seasonal changes in birth rates: when 
each species has a specific growing season disjoint from all others, and all 
resources are conservative (meaning it is a linear function of population 
densities) , a system with four resources can support more than four species 
without any one going into extinction. Our objective is to extend their ob- 
servation to competition models based on interacting particle systems. 

Recall that when two contact processes with equal death rates 5 and 
equal dispersal kernels compete with one another on Z 2 , one of the two 
species outcompetes the other one. To search for coexistence, we follow the 
work of Armstrong and McGehee and consider an environment with seasonal 
changes in birth rates: time will be divided into seasons of length D, called 
alternately season 1 and season 2. For simplicity and because this is the case 
of interest in spatial ecology, we assume that d = 2. However, all our results 
extend easily to any dimension d > 1. Since we will deal with long range 
interactions, it is convenient to assume that particles evolve on the rescaled 
lattice Z 2 /L where I is a large integer which is referred to as the range 
of the interactions. Our model is a continuous-time Markov process whose 
state at time t is a function that maps the rescaled lattice into {0, 1,2}. 
To formulate the dynamics, we let 

M{x) = {y G I?/L :0<\\x- y\\oo < 1} 

be the interaction neighborhood of site x 6 Z 2 /L, and 

i? / a\ \{y£Af(x):£(y)=i}\ 

M \Af(x)\ 

be the fraction of neighbors of site x in state i. The state of site x £ Z 2 /L 
flips according to the following transition rates: 

transition season 1 season 2 

0-1 0ufi(x,O Pufi(x,0 

0- 2 /3 2 i/ 2 (z,0 /W 2 (x,0 

1- 6i 5 1 

2- 6 2 6 2 

This indicates that each species behaves individually like a contact process. 
Their associated birth rates, however, depend on the season, making our 
process a generalization of the Neuhauser's competing model in which the 
birth rates are functions of time. We call the process the periodic compe- 
tition model. We assume for simplicity that each of the species has the same 
death rate in both seasons. However, the reduction to equal death rates is 
not important for the proofs. In particular, Theorem 1 and Corollary 3 be- 
low can be extended to the case when the death rates may vary provided 
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one replaces the death rate by the average of the death rates in season 1 and 
season 2, respectively. We now give a sufficient condition for coexistence. 

To begin with, we observe that when the dispersal range L tends to infin- 
ity, the population dynamics (in a finite volume) converges to the following 
mean-field model: 

= ^(tfutil -ui- u 2 ) - 8\u\, 

(2) ^ 

— = (3 2 (t)u 2 (l - u\- u 2 ) - S 2 u 2 , 

where Ui, i = 1,2, denotes the population density of species i and /3 % {t) are 
periodic step functions with period 2D given by 

Pn, whenO<t<A 
.$2, when D < t < 2D. 

In the absence of species 1, species 2 evolves according to the pair of logistic 
equations 



(3) 



du2 _ f 02iu 2 (l - u 2 ) - 5 2 u 2 , if t G [2nD, (2n + 1)D), 
dt I (3 22 u 2 (l - u 2 ) - 5 2 u 2 , if t £ [(2n + 1)1), (2 ? i + 2)D). 



Denoting by U2(t) the solution, the function t >— > u 2 (2nD + 1) for < t < 2D 
converges uniformly to a piecewise-smooth periodic curve U2- This follows 
from the monotonicity of the solution with respect to the initial condition. 
Moreover, since u 2 is periodic and continuous, standard topological argu- 
ments imply that the convergence is uniform on the real line. We call u 2 the 
equilibrium curve for the 2's. One can define u\ for the l's in a similar way. 
The top panel of Figure 1 gives an illustration of the density of 2's in the 
absence of l's for the interacting particle system which, in view of the size 
of the universe and the range of the interactions, is a good approximation 
of u 2 . 

We now return to the spatial model. When the birth rates (5ij are chosen 
appropriately, one can compare the l's (resp., the 2's) with a supercritical 
branching process and show that the l's and 2's coexist. Based on Armstrong 
and McGehee (1976), we expect that any number of species can coexist 
provided there are as many seasons as there are species, but for simplicity 
we restrict ourselves to a system with only two species. 

Theorem 1. Suppose that 

1 r 2D 1 r 2D 

— P\t)(l-u 2 (t))dt>5 1 and — /? 2 (t)(l-ui(t))dt>5 2 . 
ID Jo ZD Jo 

Then, when L is large, coexistence occurs. 
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density of species 2 {fi%\ = 1 and $22 — 3) 














density of species 1 (Pu — 3 and 0i2 — I) 
density of species 3 (/?3i — 2 and (3w = 2) 



density of species 2 (/?2i — 1 and P22 — 3) 




Fig. 1. Evolution of the densities of species for the interacting particle system along 4 
periods = 8 seasons. The universe is the 400 x 400 lattice with periodic boundary condi- 
tions. The length of each season is D — 10, the range of the interactions L = 200 and the 
death rates Si — 1. Our pictures show species 2 alone, species 1 and 2, and species 1, 2 and 
3 in competition, respectively. 



In the previous theorem and after, coexistence for the interacting particle 
system means that starting from any initial configuration with infinitely 
many l's and infinitely many 2's, 

liminf P(&(a?) = i) > for i = 1, 2 and all x G 1?/L. 

Note that due to temporal heterogeneities there is no stationary distribution, 
but there is an initial distribution so that the states at times and 2D have 
the same distribution. It seems reasonable to call this a periodic distribution. 
Since J(l — U2)/2Ddt is the average number of empty sites available for 
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invasion, the first integral in Theorem 1 represents the growth rate of l's 
when rate averaged over time. Therefore, if the average growth rate is greater 
than the death rate then the l's can invade a community of 2's. Similarly, 
the second inequality says that if the average growth rate for species 2 is 
greater than the death rate then the 2's can invade a community of l's. 

As previously mentioned, it has already been proved by Armstrong and 
McGehee (1976) that temporal heterogeneity promotes coexistence. Their 
result, however, was based on a deterministic nonspatial model. It is im- 
portant to extend their result to a model including both temporal hetero- 
geneity and spatial structure in the form of local interactions because the 
spatial component is identified in ecology as an important factor in how 
communities are shaped. Even though the sufficient conditions we found for 
coexistence of two species in the presence of two seasons turn out to be 
the same for both the interacting particle system and its nonspatial ana- 
log, our main result is interesting for ecologists since it is known from past 
research that spatial models can result in predictions that differ from non- 
spatial models. For instance, it has been proved theoretically in Neuhauser 
and Pacala (1999) and Lanchier and Neuhauser (2009) that the inclusion of 
local interactions may translate into a reduction of the coexistence region. 
Their models were based on particle systems with short range interactions. 
Interestingly, Neuhauser (1994) also proved that a sexual reproduction pro- 
cess with a long range dispersal kernel does not exhibit the same features 
as the corresponding mean- field model, namely, for intermediate values of 
the birth rate survival occurs for the mean-field model while the interacting 
particle system exhibits a metastable behavior and eventually goes extinct. 
In other words, even in the presence of long range interactions, spatial and 
nonspatial models may have significantly different behaviors. 

The second panel of Figure 1 shows the evolution of the densities of species 
1 and 2 in competition, which illustrates Theorem 1. The bottom panel sug- 
gests that temporal heterogeneity can also promote coexistence of three 
species even in an environment with only two seasons. Interestingly, numer- 
ical simulations also suggest that, with the birth rates indicated in Figure 
1, the three species coexist whereas in the absence of species 2, species 3 
outcompetes species 1. This indicates that the inclusion of an additional 
species may promote coexistence. In conclusion, we state the following. 

Conjecture 2. In the presence of two seasons, three species can coex- 
ist. 

The conjecture is probably very difficult to prove. In principle it can 
be approached with invadability ideas as in Durrett (2002). Suppose for 
simplicity that we have found three species that can coexist in pairs. To show 
that all three can coexist, we need to (i) show that the two species systems 
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converge to equilibrium, and that (ii) in all three cases the third species 
can invade the other two in their equilibrium. Problem (i) seems difficult. 
Indeed, we do not see how to prove that the mean field ODE converges to 
equilibrium. Due to the time inhomogeneity, one does not have Lyapunov 
functions. 

We now derive from the two integrals in Theorem 1 a (stronger but) 
more explicit sufficient condition for the coexistence of both species in the 
spatially explicit model. This condition is given in Corollary 3 below. First 
of all, we set 

P2j = lim U2(2nD + (j — 1)D) while holding u± = 

J n—>oo 

for j = 1,2. That is, p2j is the limiting density of the 2's at the beginning 
of season j in the absence of l's. Since the 2's grow according to the pair of 
logistic equations (3), the densities P21 and P22 are related by 

p(p2i,D;P2i,5 2 ) =P22 and p(p 2 2, D; (3 2 2, 82) = P21, 

where the function p is defined by 

( + a x\ n exp(rt) 

p(uo,t;P,6) = - — - — —-, 

1 — K i no(l — exp(rij) 

where 

r = /3-5 and K = l-6//3 

are the intrinsic rate of growth and the carrying capacity, respectively. By 
letting j = (nmod2) + 1, one can define U2 recursively, by 

U2(0)=P2i and u 2 (nD + 6) = p(u 2 (nD),6; fcj,^) 

for < 9 < D. One can define pn, p\2 and u\ for the l's in a similar way. 
Now, assuming that season 2 is the growing season for species 2, that is, 
/?2i < /322> the limiting curve U2 is convex and decreasing in season 1 and 
concave and increasing in season 2. (See the top panel of Figure 1 for an 
approximation of U2-) In particular, 

■j- 1° l-u 2 {t)dt>l- P2l+ J 22 and 1 l-u 2 {t)dt> I-P21. 
V Jo 2 D Jd 

This indicates that the densities of empty sites available for the l's in seasons 
1 and 2 are bounded from below by 1 — (p2i +P22)/2 and 1— P21, respectively. 
Therefore, the l's dominate a contact process with birth rate 

0u(l- P -^^) and M1-P21) 
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in seasons 1 and 2, respectively. Similarly, assuming that 0n > 0i2, the 2's 
dominate a contact process with birth rate 

Pn+Pn 



02l(l-?>l2) and /?22^1 

in seasons 1 and 2, respectively. The symmetry is flipped because the 2's 
have a different growing season. By applying Theorem 1, we finally obtain 

Corollary 3. Suppose that 0n > 0i2 and 022 > 021, an d that 
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1 



P21 + P22 



#2l(l -P12) +&2 1 



+ {3 12 (1-P21 
P11+P12 



> fa, 



> 5 2 



where pij is the density of species i in the absence of the other species at the 
beginning of season j in the mean-field equilibrium. Then, when L is large, 
coexistence occurs. 

In the previous result, each species does better than the other in one 
season, but this is not necessary for coexistence. For a concrete example, 
suppose that 

D = l, 021 = 5.2, 022 = 1-0 and £ 21 = <5 22 = 2.0. 

Solving the equations numerically shows that the integral of the density 
U2(t) over the first season is approximately 0.366066. Thus, if we pick 

0ii = 10000, 012 = 0, $11 = 6000 and 5 12 = 100, 

then we have 0ii/<5n = 5/3<13/5 = 02i/$2i but 

1 f 2D #11 f 1 

— [3 1 (t)(l-u 2 (t))dt = ^ (1 -u 2 {t))dt> 3169 > 3050 
ID Jo 2 Jo 



Sn+S 
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so the l's can invade the 2's. In other respects, 

ui(t) < l-<5n/0ii = 0.4 fort<l, 

«i(t) < 0.4exp(-100(t - 1)) for 1 < t < 2, 

which implies that 



1 
2D 



21) 



P 2 (t)(l-u 1 (t))dt- 



0: 



21 



/ (1- Ul (t))dt + 
J 



022 

2 



(l-«i(i))di 



0: 



3 498 

io A ' + Iooo ft2 = 2058> 



exp(-100t)dt 
$21 + $22 
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so the 2's can invade the l's as well. The fact that coexistence may occur 
even though species 2 is a better competitor than species 1 in both seasons is 
due to the combination of the presence of seasons (temporal heterogeneity) 
and the fact that species 1 evolves at a much faster rate than species 2. In 
our example, the density of l's immediately goes to the carrying capacity in 
season 1 and then drop quickly to in season 2. Individuals of species 2, on 
the contrary, evolve slowly so they never really get to full strength: Their 
effective average density is less than in the homogeneous case, which allows 
the l's to invade the 2's, even thought the 2's are superior competitors. 

2. Mean-field model. In this section, we prove a coexistence result for 
the mean-field model which is analogous to Theorem 1. A common way for 
showing coexistence of two species models is to prove that both species are 
mutually invadable, that is, species 1 can invade a community of 2's when the 
l's are small in number, and vice versa. This implies that their densities are 
bounded away from the axes on the u\-u 2 plane. In Proposition 2.1, below, 
coexistence for the mean-field model means that starting with a positive 
density of type 1 and type 2, 

liminf Ui{t) > fori = 1,2. 

Note that, similarly to the interacting particle system, due to temporal het- 
erogeneities, the mean-field model (2) has no fixed point but instead con- 
verges to a limit cycle. 

Proposition 2.1 (Coexistence). Suppose that 

l r 2D l r 2D 

— /3 1 (t)(l-u 2 {t))dt>S 1 and — (3 2 (t)(l - u^t)) dt > 6 2 . 

ID Jo ZD Jo 

Then there is coexistence in the mean-field model (2). 

Proof. By symmetry, we only need to prove survival of the l's. The 
idea is to show the existence of a small e > 0, such that if U\(t) € (0,e) for 
all t G [0, 2D] and 

/ t3 1 (t)(l-u 2 (t))-6 1 dt>0, 
Jo 

then u\{2D) > ■ui(O). Assuming that ui(t) £ (0,e), we obtain 

^ >^ 1 (t)ni(l -e-u 2 ) - Sim. 
Dividing both sides by u\ and integrating from to 2D gives 

log«i(2D)-log« 1 (0)>(l-e)(/3n + / 9i 2 )D- / (3 1 (t)u 2 (t) dt - 26^ . 

Jo 
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Therefore, if u 2 satisfies the condition 

(4) (fin +Pw)D- p 1 (t)u 2 (t)dt>25 1 D 



o 

one can choose e > small so that u\(2D) > ui(0). If 1x2(0) < P21 then 
v>2(t) < u 2 (t) at all times, which implies that 

(Pii+Pi2)D- PHt)u 2 {t)dt> {3 1 {t){l-u 2 {t))dt>25 l D, 
Jo Jo 

so (4) is satisfied. This proves the result when 112(0) < p 2 ±. If on the contrary 
1*2(0) > p 2 i, then, by monotonicity of the solution with respect to the initial 
condition, survival of type 2 holds as well. This completes the proof. □ 

Note that taking u\(t) = in Proposition 2.1 gives the following sufficient 
condition for the survival of species 2 in the absence of species 1. 

Corollary 4 (Survival of a single species). In the absence of 1's, 

liminfii2(i) > whenever @ 22 > 

t^oo 2 

3. Preliminaries on periodic contact processes. From now on, we call 

a periodic contact process a contact process on Z 2 /L in which particles 
give birth at rate a(t) and die at rate 5, where a(t) is positive, piecewise 
continuous and periodic. A particle produced at site x £ Z 2 /L is sent to 
a site chosen uniformly from J\f(x). The birth occurs if the target site is 
empty, and is suppressed otherwise. This process is denoted by A t (a,5) and 
thought of as a set valued process, that is, A t (a,5) is the subset of 1?/L 
corresponding to the set of sites which are occupied at time t. The dual 
process of the periodic contact process starting at site x at time t is denoted 
by Ag' 1 and represents the descendants of (x,t) at time t — s. Note that 
contrary to the basic contact process, the periodic contact process is not 
self-dual since particles in the dual process A^'* give birth at rate a(t — s), 
which a priori is different from a(s). 

Similarly, we call a periodic branching random walk on M 2 a branching 
random walk in which particles give birth at rate a(t) and die at rate 5. 
A particle produced at x G M 2 is sent to a site chosen uniformly at random 
from x+ [—1, l] 2 . This process is denoted by Zt(a,S). For both the periodic 
contact process and the periodic branching random walk, we use a super- 
script to indicate the initial configuration. We first prove some convergence 
properties of A t (a,5). 

Lemma 3.1 (Convergence to branching random walk). Let A = {0} and 
T > 0. Then, as the range L^oo, the periodic contact process Af(a,5) 
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converges weakly to Z^(a, S) in the Skorohod space of cadlag functions from 
[0,T] to the subsets ofR 2 . 

Proof. This is a classical result when the birth rate a(t) is a constant 
[see, e.g., Durrett (1991) for a proof]. Let b = maxa(t). Then, for any t <T, 
the number of births by time t in the contact process is dominated by the 
number of particles in Z^(b,0) from which it follows that 

E\Af{a,5)\ <E\Zf(b,0)\ =e bt <e bT . 

Then, by Markov's inequality, we have 

P{\A?{a, 5)\ > L 1 / 3 ) < - as L ^ oo. 

Now, conditional on E\ = {\Af(a,S)\ < L 1 / 3 for all t < T}, at each birth, 
the probability that the offspring is sent to a site already occupied is less 
than 

£l/3 

LV3.___^o asL^oo 

since the size of the interaction neighborhood is bounded by (2L + l) 2 . We 
denote by Ei the event that, up to time T, all the offspring of the periodic 
contact process are sent to empty sites. Then, on the event E\ n E%, one can 
create a coupling in which all the particles in Af(a,8) are within distance 
L x / 3 /L in the uniform norm of their counterparts in the branching random 
walk. For both processes, we call particle the particle at site at time 0, 
and particle i, i > 1, the particle produced at the ith birth event. For i£l 2 , 
let ttl(x) be the (unique) vector such that 

\\x-ttl(x)\\ 00 <1/2L and tt l (x) eZ 2 fL. 

The processes are coupled as follows. If particle i in the periodic branching 
random walk dies at time t$, then particle i in the periodic contact process 
dies at time ti as well. Now, assume that particle i in the periodic branching 
random walk is located at x and that particle % in the periodic contact 
process is located at xl- If particle i in the periodic branching random walk 
produces an offspring which is sent to y at time tj, then particle i in the 
periodic contact process produces an offspring sent to xl + ^l{v — %) at time 
ti. Since there are at most L 1 / 3 birth events by time T, adding the errors we 
find that, on the event EiHE^, the distance between particles in the contact 
process and their counterparts in the branching random walk is bounded by 

L 1/3 sup vr L (x) = L 1/3 /2L — ► as L — > oo. 
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Since P{E\ f] E 2 ) — ► 1 as L — > 00, we conclude that the periodic contact 
process Af(a,5) converges weakly to the periodic branching random walk 
Z^-(a, 6) as L — ► 00. □ 

We now prove that the number of particles in the periodic contact process 
At(a,5) is almost deterministic when L is large: If we focus our attention to 
a space-time box of finite size, one can show that the number of particles 
does not deviate much from its expected value. 

Lemma 3.2 (Convergence to expected value). Assume that a(t) is piece- 
wise constant and periodic. For AcM 2 , let Af denote the periodic contact 
process starting from An (Z 2 /L) . For any subset A, any site x on the rescaled 
lattice 1? /L and any time t > 0, let 

N(x,t) = \A?C\{x+[Q,l] 2 }\. 

FixT>0 and S>0. Then for all e>0, 

P(\N(x,t)-EN(x,t)\ >5L 2 e for some 

x G [—Vt, Vt] 2 and some t < S) — > as L —> 00. 

Proof. We follow the proof of Lemma 3.5 in Durrett and Lanchier 
(2008), in which the result is proved when the birth rate a(t) is a constant. 
Let 

N £ (x,t) = \A t i n{x + [0,e] 2 }\ and y,z £ x + [0,e] 2 D (Z 2 /L). 
First, we observe that the number of descendants in the dual processes 

^ Z t 

and A t ' depends on t and the birth and death rates, but not on L. In 
particular, since the interaction neighborhood contains (2L + l) 2 — 1 sites, 
the probability that the dual processes hit, that is, that a descendant is sent 
to a site already occupied by another descendant, is < C\/L 2 . Now, noting 
that 

. A . \ fo, if Af l nA = 0, , . A . , fo, it Ai ,t nA = 0, 

Af(y) = I ' * 'and Af(z) = < * 

tV> \l, ifAfnA^0 tK ' \l, ifAfnA^0, 

a standard construction shows that the covariance of Af-(y) and A^(z) can 
be bounded by the probability that the dual processes hit [see page 21 in 
Griffeath (1979)]: 

(5) co V (At(y),At(z))<C 1 /L 2 . 
Since N £ (x,t) = J2y£x+[o,e] 2 Af(y), inequality (5) implies that 
Y&rN £ (x,t) < e 2 L 2 + {e 2 L 2 ) 2 Ci/L 2 < C 2 L 2 e 2 . 
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Then, by Chebyshev's inequality, we obtain 

(6) P{\N £ (x,t)-EN £ (x,t)\>L 2 e 3 )<C 2 /e 4 L 2 ^0 asL^oo. 

This holds for all x € [~VT, VT + l] 2 n (Z 2 /L) and all t < S. To deduce 
the result, we first use a discretization of space and time. The idea is to 
rely on the convergence in (6) to control the densities N £ (x,t) at a finite 
number (that does not depend on L) of space-time points. More precisely, 
we let t > be a positive constant and look at the densities N £ (x,t) for 
space-time points (x,t) belonging to a subset of eZ 2 x rZ + . As L — ► oo, 

P(\N £ {x,nT) - EN £ (x,nr)\ > L 2 e 3 for some < n < m and 

some x G [-y/T, VT + l] 2 n el?) 

<C 3 T(m+l)/e 6 L 2 0, 

where m = min{n > 1 : nr > S 1 }. Using large deviation estimates for the Bi- 
nomial distribution and taking r > small, one can show that the number 
of sites in x + [0, e] 2 that flip between nr and (n + l)r is smaller than 2L 2 e 3 
with probability at least 1 — Csexp(— aL 2 ) so that 

P(\N £ (x,t) - EN £ (x,t)\ > 3L 2 e 3 for some t < S 

and some x E [—VT, y/T + l] 2 n eZ 2 ) — > as L — >oo. 

To conclude, we observe that if 

|iV e (x,t) -SiV e (x,t)| <3L 2 e 3 for all x € [- \/f, + l] 2 n eZ 2 , 

then we have that 

\N(x, t) - EN{x, t) | < 3L 2 e 3 /e 2 + 2L 2 e = 5L 2 e 

for all x € [— VT, VT] 2 n (Z 2 /L) since there are at most e -2 squares with 
length side e included in the unit square x + [0, l] 2 and at most 2L 2 e re- 
maining sites. □ 

Lemma 3.3 (Convergence to ODE). Let A t =A t (/3 1 ,6) be the periodic 
contact process starting from the "all occupied" configuration. For any site 
x on the rescaled lattice Z 2 /L, we set 

u L (x,t) = P(xGA t ((3 1 ,5)). 

Let T > 0. Then, as the range of the interactions L — > oo, the function t ^ 
ul{x,€) converges uniformly on [0, T] to a solution u(t) of the ODE 

'^L=f3 1 (t)u{l-u)-5u withu(0) = l. 
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Proof. Since the evolution rules are translation invariant, the probabil- 
ity of {x £ A t } does not depend on x. Therefore, we can set UL(x,t) =UL(t). 
Let T > and S(t) be the semigroup that generates the process At- For 
any cylindric function /, that is any function that depend on finitely many 
coordinates, the Hille-Yosida theorem states that 

(7) j t s(t)f = smf, 

where f2 is the Markov generator of S(t). Now, letting f(A t ) = ^{xeA t }i we 
have 

(8) S (t)f = Ef(A t ) = P(xGA t ) = u L (t). 
Combining (7) and (8) implies that 

(9) ^-P(xeA t )= tj?t%P(x ^ A t and ye A t ) -SP(xe A t ). 

We claim that {x £ At} and {y £ At} are asymptotically independent for 
any t <T, that is, 

(10) P(x $ A t and y £ A t ) -> P(x A t ) ■ P{y £ A t ) as L -> oo. 
To see this, we first observe that (10) is equivalent to 

(11) P{A X / = and Af + 0) -> P(i^* = 0) • P(if * / 0) as L -» oo. 

Since t <T, the number of particles in AJ'* and A|' , s <t, depends on T and 
the birth and death rates but not on L so, with probability arbitrarily close 
to 1 when L is sufficiently large, the number of particles is bounded by L 1 / 3 
(see the proof of Lemma 3.1). Moreover, on such an event, the probability 
that the dual processes hit is less than 2L 1//3 /4L 2 — ► 0. In particular, the 
dual processes are independent in the limit as L — > oo which establishes (11) 
and (10). Finally, combining (9) and (10) gives the desired result. □ 

4. Proof of Theorem 1. 

4.1. Block construction. The first step to prove Theorem 1 is to construct 
the periodic competition model £t from a Harris' graphical representation 
[Harris (1972)]. For x,y£Z 2 /L, wc let 

(\Af(x)\-\ if 0<||x-y|| oo <l, 
\ 0, otherwise. 

For any ordered pair of sites (x,y) and any i,j £ {1,2}, we introduce the 
following collections of independent Poisson processes: 

arrival times rates symbols 

T n( x >y) Pij<l(x,y) x-^y 

Un(x) S{ x at x 



q(x,y) 
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that is, we draw an i-arrow from site x to site y at time T^(x,y) to indicate 
that if x is occupied by species i, site y is empty, and 

(12) \I%(x,y)/D]:=mm{aeZ:a>I%(x,y)/D} = j (mod 2) 

then site y becomes occupied by an individual of species i. Note that the 
condition in (12) means that the arrival time T^(x,y) occurs in season j. In 
addition, we put a x at site x at time U^(x) to indicate that an individual 
of species i at x dies. A result of Harris (1972) implies that the previous 
construction induces a well-defined Markov process. Moreover, the flip rates 
indicate that this process is the periodic competition model. 

Our proof of Theorem 1 relies on the use of a block construction, which 
is now a standard technique. The idea is to show that, when viewed on 
suitable length and time scales, our process dominates the set of wet sites 
of a Af -dependent oriented percolation process on 

£ = { (m, n) G Z 2 : m + n is even and n > 0} 

in which sites are open with probability p close to 1. More precisely, each 
site (m,n) G C is associated with a random variable u{m,n) G {0,1} which 
indicates whether the site is open (state 1) or closed (state 0). The M- 
dependency means that 

P(uj(mi,rii) = 1 for 1 < i < k) = p k , 

whenever \\(mi,rii) — (mj,rij)\\ > M for i^j. A site (m,n) is said to be wet 
if there exists a sequence of integers mo, mi, . . . , m n = m such that: 

1. For i = 0, 1, . . . ,n — 1, we have — rrii\\ = 1. 

2. For i = 0, 1, . . . ,n, site (mi, i) is open, that is w(mj, i) = 1. 

See Durrett (1995) for more details. To compare the periodic competition 
model viewed on suitable length and time scales with oriented percolation, 
we let 7 be a large parameter to be fixed later, and let Ni(x,t) denote the 
number of type i particles in x + [0, l] 2 at time t. Assuming that 

(13) Ai(0,0) > L 2 exp(- 7 T) fori = 1,2, 

we will prove that with probability close to 1 when the parameters T and L 
are large, 

(14) Ni((s/T,0),T) >L 2 exp(- 7 T) fori = 1,2. 
In particular, if we call (m,n) G C an occupied site whenever 

Ni(m(VT,0),nT) >L 2 exp(- 7 T) for i = 1,2, 

then the set of occupied sites dominates the set of wet sites of an oriented 
percolation process in which sites are open with probability p close to 1. 
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The construction of the percolation process relies on a coupling argument in 
which the random variables u(m,n) are defined by induction over the level 
n. This technique is now standard and we refer the reader to the Appendix 
of Durrett (1995) for a detailed construction of the process in the general 
case. Since percolation occurs whenever the density p is sufficiently close to 
1 (i.e., starting with one open site at level 0, there is a positive probability 
that the cluster of wet sites is infinite) , Theorem 1 follows from the compar- 
ison with oriented percolation. In conclusion, the aim is to prove that the 
conditional probability of the events in (14) given the events in (13) can be 
made arbitrarily close to 1. 

Since the proofs of the inequalities in (14) for i = 1 and i = 2 are identical, 
we will only show the first one (i = 1). There are three steps in showing 
the survival of species 1 (see Figure 2 for an illustration of the three steps 
described below): 

(i) The density of type 2 particles in [— y/T, \/T] 2 between times \/T and 
T is (almost) bounded by the mean-field equilibrium U2(t). 

(ii) If the unit square [0, l] 2 contains at least L 2 exp(— 7T) type 1 particles 
at time 0, then it will contain at least L 2 exp(— "yT + cT) type 1 particles 
at time T — 2\/T for some c> 0. 

(iii) If the unit square [0,1] 2 contains at least L? exp(— 7T + cT) type 1 
particles at T — 2\/T, then at least L 2 exp(— 7T) of them will "invade" 
the unit square (VT, 0) + [0, l] 2 at time T. 

4.2. Proofs of (i)-(in). We start with condition (i) which is proved in 
Lemma 4.1. 

Lemma 4.1. For any e > and for T sufficiently large, 
P(N 2 (x,t) > L 2 {u 2 (t) + 7e) for some 

x G [—y/T, VT] 2 and some t G [VT,T]) -^0 as L — > 00. 

Proof. Let At be the periodic contact process with birth rate /3 2 (t), 
death rate 62, and initial configuration Aq = 2, and let £t = {x:£t(x) = 2} 
denote the set of sites occupied by a type 2 particle at time t in the periodic 
competition model. Then A t and £t can be coupled in such a way that Q C A t 
for all t > 0. In particular, it suffices to prove the result for 

N(x,t) = \A t D{x + [0,l} 2 }\ 

instead of N2(x,t). Let u(t) be the solution of the ODE 



— = 2 {t)u{l - u) - 5 2 u withu(0) = l. 
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Fig. 2. Boxes in the block construction for Theorem 1. 

Then there exists T sufficiently large such that 

(15) |u(*)-u 2 (t)|<e foraUt>Vf. 

Moreover, by Lemma 3.3, P(x G At) — ► uniformly on [0, T] as L — ► oo 
which, together with inequality (15), implies that 

EN(x,t) <L 2 (u 2 (t) + 2e) 

(16) 

for all G [--n/T, VT] 2 x [Vt,T] as L -> oo. 
In other respects, by Lemma 3.2, 

P(|JV(z,t) -£7JV(x,t)| > 5L 2 e for some 

(17) 

x G [— VT, VT] and some t G [VT, T]) — ► as L — > oo. 
The result follows by combining (16) and (17). □ 

Let Zt = Zt(a,S) be the restriction of the branching random walk with 
birth rate a(t) and death rate 5 to the square It = [— VT, VT] 2 , that is, 
particles landed outside It are killed. In preparation for proving (ii), we 
need to show that if the averaged birth rate a(t) along one period is greater 
than the death rate 5 then the number of particles in Z% grows exponentially 
fast. 
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Lemma 4.2. Assume that Z = Z§ = {x}. For s > r > fixed and T 
sufficiently large 

E\Zf n {y + [0, 1] 2 }| > exp Q j\a{6) - 5) d0 
for all x, ye [- VT/2, VT/2} 2 and t G [rT,sT]. 

Proof. Let Sf = Sf{a) be the random walk that jumps at rate a(t) 
and has kernel the uniform distribution on the square [—1, l] 2 , that is, the 
new position of the process is chosen uniformly at random from Sf + [— 1, l] 2 , 
and denote by Sf its restriction to the square It- Let 

m(t,x,A)=E\Zf(a,5)f]A\ 

denote the mean number of particles in the set A at time t. Then we claim 
that 

(18) m(t,x,A) =exp(j\a(0) - 5)d0)p(Sf(a) e A). 

This holds because both sides of (18) satisfy the differential equation 
dm(t,x,A) f 

— = —dm(t,x,A) + / a(t)m(t ) x,dy)l^_ y / Tiy /j> i ayi , {A-y), 

where v is the uniform probability measure on [—1, l] 2 . Since < mma(t) < 
maxa(t) < oo, one can rescale time so that Sf has constant jump rate 
min a(t). Let 



t = uF(T), x = v^/F(T), y = w^F(T), v,w £ [-i, |] 2 . 

Here F(T) = 0(T) is piecewise continuous with 1 < F{T)/T < maxa(t). An 
application of the invariance principle shows that the process Sf converges in 
probability to a Brownian motion when time and space are properly scaled. 
More precisely, sending T — > oo we have 

F(T)-P(Sf£y+[0,l] 2 )^p u (v,w), 

where p u (v,w) is the time-n transition probability of a 2-dimensional Brow- 
nian motion from v to w that is killed outside [— 1,1] 2 . Since \v — w\ < 1, 
the probability p u (v,w) is bounded from below which implies that, for T 
sufficiently large, 

P(Sf ey+[0,l] 2 ) = O(T- 1 ). 

It follows that 

m(t, x, y + [0, l] 2 ) = exp (J (a{6) - 5) d9^j ^(T^expQ J* (a(d) - 5) d6 
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for all t G [rT, sT] with T sufficiently large. □ 

In the next lemma, we establish (ii). The idea is that, since the density 
of type 2 particles is bounded by ui + 7e, the l's grow exponentially fast 
from time VT to time T — 2y/T so that they can regain all possible losses 
occurring by time VT. 

Lemma 4.3. Assume that iVi(0,0) > L 2 exp(-7T). Then, for 7 and T 
sufficiently large, 

P(Ni(0, T - 2VT) < L 2 exp(-7T + cT)) as L -> cx> 

for a suitable constant c > 0. 

Proof. First of all, we observe that the probability that a type 1 particle 
remains alive VT units of time is given by p = exp(— <5i\/r). Therefore, an 
application of the large deviation inequality for the Binomial distribution 
implies that as L — > 00, 

P(iVi(0, VT) < {L 2 /2) exp(- 7 T - SiVf)) < exp{-e~~< T ~ SlVT L 2 /8) 0. 

In particular, since the inclusion relation C on the set of sites occupied by 
type 1 particles is preserved by the dynamics, it suffices to prove the result 
when 

N 1 (0,VT) = (L 2 /2)exp(-~{T-5 1 VT) and 

(19) 

^(x)^=l for all xi [0,1] 2 . 
The assumptions of Theorem 1 allow to fix e > such that 

(20) / /3 1 (t)(l-u 2 (t)-8e)-5i(it>0. 
Jo 

Applying (18) to the restricted branching random walk Z-t with birth rate 
a(t) = (3 l {t) and death rate 5 = 5\, which dominates the number of type 1 
particles in the square It, and taking (19) into account, we obtain that, for 
all sites x G It and all times t G [VT, T] , 

Mi (x,t)< iVi (0, VT) exp [J^-il 31 dO^j 

< L 2 exp(-7T + max( / 9n, / 9i2)r). 
In particular, there exists a sufficiently large 7, fixed from now on, such that 

(21) P(N 1 (x,t) < L 2 e for all x G I T and all t G [VT, T — 2VT\) -)• 1 

as L — > 00. 
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By Lemma 4.1, the density of type 2 particles inside It = [— VT, VT] 2 is 
bounded by u 2 + 7e between time y/T and time T — 2y/T. This, together 
with (21), implies that, with probability arbitrarily close to 1 when L is 
large, the number of empty sites in each square x + [0,1] 2 is larger than 
L 2 (l-u 2 (t) -8e). So, in a given interval, as L — > oo, the number of l's 
dominates the number of particles in a branching random walk with birth 
rate 

(22) a(t)=/3 1 (t)(l-u 2 {t)-8e) 

and death rate §i, and with (L 2 /2) exp(— 7 T — SiyT) particles in [0,1] 2 at 
time y/T. To see this, note first that (22) is a lower bound of the rate of 
successful birth of type 1, that is, the rate at which type 1 particles are 
actually created, in the periodic contact process when (21) holds. Moreover, 
the arguments of Lemma 3.1 imply that both processes can be coupled in 
such a way that type 1 particles in the periodic contact process are arbitrarily 
close to their counterparts in the branching random walk. In particular, it 
suffices to prove the result for the branching random walk described above 
restricted to the square It- We denote this process by Z t (a,5i) and let 

N{x,t) = \Z t (a,5 1 )n{x + [0,l} 2 }\. 

By Lemma 4.2 and (20), for T sufficiently large, 

EN(0, T - 2x/T) > (L 2 /2) exp(- 7 T - 5iVT) 

xexp(-^_ (3\t)(l-u 2 ( y t)-8e)-5 1 dtj 

> L 2 exp(- 7 T + 2cT) 

for a suitable constant c > 0. Finally, by Lemma 3.2 (which applies to the 
periodic branching random walk as well), 

N(0, T - 2VT) > L 2 exp(- 7 T + 2cT) - 5L 2 e > L 2 exp(- 7 T + cT) 

with probability arbitrarily close to 1 for L and T sufficiently large. □ 

To complete the proof of Theorem 1, it remains to "move" our parti- 
cles from the square [0, l] 2 to the square (VT,0) + [0,1] 2 and show that 
L 2 exp(— 7 T) of them will live until time T. Our last lemma shows (iii) in 
the list. 

Lemma 4.4. Assume that N± (0, T - 2y/T) > L 2 exp(- 7 T + cT) and that 
N 2 (x,t)<L 2 (u 2 (t) + 7e) for all x £[-x/T,VT] 2 and t G [T - 2\/T,T}. 
Then, for T large, P(Ni((\/T,0),T) < L 2 exp(- 7 T)) ^0 as L -> 00. 
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Proof. Fix a > 0. For any integer m£Z and any time t > 0, we let 

B m = (m/2, 0) + [0, 1/2] 2 and H m , t = |{iG B m : &(s) = 1}|. 

Assume first that H m+ ij > 2L 2 exp(— aT) for some i G (s, s + 1). The prob- 
ability that a type 1 particle remains alive at least one unit of time is e~ Sl . 
Since these events are independent, the random variable H m+ i iS+ i domi- 
nates a Binomial random variable with parameters 

n > 2L 2 exp(— aT) and p = e~ Sl . 

An application of the large deviation result for the Binomial 

P(X < n(p - z)) < exp(-nz 2 /2p) 

with z = p/2 then gives 

P(H m+1)S+1 < L 2 exp(-aT - 5i)\H m+u > 2L 2 exp(-aT) 

(23) for some t£ (s,s + l)) 

< exp(-np/8) < exp(- e - aT -' 5l L 2 /4). 

Now, assume on the contrary that H m+ ij < 21? exp(— aT) for all t £ (s, s + 
1). Then the probability that a type 1 particle in B m at time s gives birth 
by time s + 1 to an offspring which is sent to an empty site in B m+ \ and 
that both the parent and the offspring live until time s + 1 is bounded from 
below by 

Po = Jee~ 25l {l ~ e" b )(l - max(p 21 ,p 22 ) - 7e - 2exp(-oT)), 

where b = min(/?n, ^12 )• The terms on the right-hand side are lower bounds 
for the probabilities that (i) the offspring is sent to B m+ i, (ii) particles live 
until time s + 1, (iii) the birth occurs by time s + 1, and (iv) the landing 
site is empty (since U2(t) < max(p2i,P22))- This implies that 

P(H m+ i s+1 < (p /2)H m s \H m+ i t < 2L 2 exp(-aT) 

(24) 

for all t E (s, s + 1) and H m , s > L 2 exp(-jT)) -> 
as L — > 00. Applying (23) and (24) with m = 2%/?" — n and s = T — n gives 

H 2VT-n+l,T-n+l > ^ min (P0, e ~ Sl ) H 2VT-n,T-n = ^ H 2VT-n,T~n 

with probability close to 1 when L is large on the event that H mtS > L 2 exp(— 7T) . 
Assuming that one fourth of the particles of Ni(0,T — 2y/T) are found in 
H Q T _ 2y /r and observing that 

\ ( y) ^L 2 exp(- 7 T + cT) > L 2 exp(- 7 T) 
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for T large enough, a simple induction allows us to conclude that 

n 2VT,T ^ 2 2VT-1,T-1 \ 2 J 0,T-2VT 



> Z (j) L exp(—jT + cT) 

with probability close to 1 when L is large. Recalling the definition of H m j, 
we obtain that, for T sufficiently large, and with probability arbitrarily close 
to 1 when L is large, 

iVi((VT,0),r)>L 2 exp(-7T). 
This completes the proof. □ 
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